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A possible mechanism for the removal of the orbital degeneracy in RTi03 (where R=La, Y, ...) 
is considered. The calculation is based on the Kugel-Khomskii Hamiltonian for electrons residing in 
the ti g orbitals of the Ti ions, and uses a self-consistent perturbation expansion in the interaction 
between the orbital and the spin degrees of freedom. The latter are assumed to be ordered in a 
Neel state, brought about by delicate interactions that are not included in the Kugel-Khomskii 
Hamiltonian. Within our model calculations, each of the t2 g bands is found to acquire a finite, 
temperature-dependent dispersion, that lifts the orbital degeneracy. The orbital excitations are 
found to be heavily damped over a rather wide band. Consequently, they do not participate as a 
separate branch of excitations in the low-temperature thermodynamics. 

PACS numbers: 75.30.Et, 71.27.+a, 75.50.Ee. 



I. INTRODUCTION 

The Mott insulator lanthanum titanate (LaTiOs) is 
one of the most puzzling materials in the family of tran- 
sition metals (TM) of the perovskite structure (see, e.g., 
Ref. 1). Each of the Ti ions in this material has a single 
valent electron in one of the three-fold degenerate t^g or- 
bitals, causing the ground state to be enormously degen- 
erate. Whereas in the majority of the TM perovskites 
the conventional Jahn- Teller distortion lifts the orbital 
degeneracy, such distortions are believed to be small in 
the lanthanum titanate, 2 and will be assumed to be zero 
for the purposes of the present paper. Here we pose the 
following question: is there an alternative mechanism to 
the Jahn- Teller one, which is capable of removing the 
orbital degeneracy? 

This question has already come up in materials with 
doubly-degenerate e^g valent electrons. In these materi- 
als, there exists an orbitally ordered phase, with coher- 
ent excitations ('orbitons'). 1,3 These excitations were re- 
cently observed in LaMn03. 4 The notion of an orbital dis- 
ordered, liquid, state was developed in Ref. 5 for the hole 
doped LaMnC>3 with a high hole concentration, where 
the antiferromagnetic (AF) insulating state gives way to 
a ferromagnetic metallic phase. Theoretically, the or- 
bitons in doped materials are constructed similarly to 
spinons, which appear within the framework of the t — J 
model for the hole-doped cuprates. 6 Another scenario of 
quantum liquid formation was proposed for the insulat- 
ing KCUF3. 7 According to Ref. 7, the orbital fluctua- 
tions (random changes of the bond directions) may be 
the source for the "quantum melting" of the magnetic 
order, because of the strong interaction between orbital 
and spin degrees of freedom. Melting, in this scenario, 
means the disappearance of a long-range AF order in fa- 
vor of a spin liquid state. 

The lifting of the degeneracy of the t^g orbital triplet 



in LaTiC>3 is a more difficult problem. In this 
cording to the analysis of Kugel and Khomskii (KK), 8 
the spin and the orbital degrees of freedom are not sep- 
arated. Hence, the spin excitations are necessarily in- 
timately involved in the formation of the orbital exci- 
tations and vice versa. A possible scenario of orbital 
liquid formation in LaTiC>3 has been proposed in Ref. 
9. These authors have conjectured the appearance of a 
resonance between the spin and the orbital excitations. 
To describe this resonance, they have used mean-field de- 
coupling, generalizing the approach proposed in Ref. 6 for 
the spin liquid state in the two-dimensional t — J model. 
This decoupling leads inevitably to the conclusion that 
the resonance states should give a contribution linear in 
the temperature to the low-energy specific heat. Such 
a Sommcrfeld-like contribution signals the lifting of the 
orbital degeneracy by forming a quantum liquid. Unfor- 
tunately, a direct experimental check of this prediction 
in LaTiC>3 has yielded negative results. 10 

In the present paper, the main idea of Khaliullin and 
Maekawa 9 is rehabilitated by means of a more refined 
approach, which goes beyond the mean-field decoupling. 
We use the methods recently tested in the theory of 
Kondo lattices, where the spin liquid state is described 
as short-range correlations of the RVB type, which arise 
due to a non-local exchange through the Kondo screen- 
ing clouds between the magnetic ions. 11 We show that 
the two-magnon excitations play a similar role in the for- 
mation of the orbital liquid correlations. 

Our analysis is based on the KK Hamiltonian 8 (see 
below). In its simplest symmetric form, this Hamilto- 
nian does not support long-range magnetic order at any 
non-zero temperature, due to some hidden symmetries. 12 
However, the Hamiltonian of the real LaTiC>3 contains 
many additional (possibly small) terms, which restore 
long range AF spin order below some Neel temperature 
Tn, as indeed observed experimentally. 2 In the follow- 
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ing, we therefore assume the existence of this ordered 
magnetic state, and concentrate on its effects on the low 
temperature dynamics of the orbital degrees of freedom. 
Once the magnetic order is stabilized, we expect these 
dynamics at temperatures T -C T/v to be described by 
the KK Hamiltonian, as discussed below. Indeed, we find 
heavily damped orbital excitations, which form a contin- 
uum whose energy scale is the same as that of the spin 
excitations, such that both orbital and spin excitations 
merge into a common continuum. As a result, a direct 
experimental observation of the orbital dynamics alone, 
at low temperatures, is hardly possible. 

The paper is organized as follows. In the next section 
we review the derivation of the KK Hamiltonian, and set 
the basis for our perturbation expansion. The latter is 
carried out in Sec. Ill, where the orbital spectrum is 
obtained and discussed. Section IV includes our conclu- 
sions. 



i 

where is the spin operator on site i, and 
J = 2t 2 /U 



(4) 
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sets the energy scale of the spin excitations. The orbital 
degrees of freedom are now described by the operators 
a\ (b\, cj), which create a spinless fermion (termed 'or- 
biton') in the yz (zx, xy) orbital on site i. The intermin- 
gling of the orbiton degrees of freedom with those of the 
spin is contained in the second term of Eq. (4) , in which 
KK introduced the orbiton operator product Y z , given 
by 



alaiOjCij 



(6) 



II. THE KUGEL-KHOMSKII HAMILTONIAN 

Following the KK approach, 8 we consider a simplified 
version of the titanate, concentrating only on the sublat- 
tice of Ti ions. Each Ti ion possesses a single electron 
in one of the three-fold degenerate t^g orbitals of the d- 
shell. Including only a single Coulomb-Hubbard on-site 
repulsion on the Ti ion, with energy U, the simple initial 
Hamiltonian reads 



H 



U 



^ ^ ^ ^ ^ ^ ^ima^im'a' (1 ^mm'^crcr') 
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where (ij) denotes nearest-neighbor pairs. Here, d\ ma 
creates an electron of spin a at site i, in one of the 
i 2 g levels, denoted m, and n ima is the number opera- 
tor in the m orbital on site i with spin index a. Kugcl 
and Khomskii 8 treated this Hamiltonian in second-order 
perturbation theory, assuming that the overlap integrals 
tYi m are much smaller than U , in order to obtain an effec- 
tive Hamiltonian describing the spin interactions among 
the Ti ions. Operating within strictly cubic symmetry, 
KK assumed that 



The expressions for TL^ K (for bonds along the re- 
direction) and for H^ K (for bonds along the y-direction) 
are derived from Eq. (4) by cyclic permutations of the 
operators a, b, and c. The full exchange Hamiltonian is 
the sum of the three, 
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Then, the first term in Eq. (4), when added together 
with the corresponding terms in H^ K and yields 
twice the total number of electrons in the system, 



n = Y,( a l a * +b l b * + c l c i)> 



(8) 



which is just a constant. 

Contrary to the situation in the manganites, the spin 
and charge degrees of freedom in the titanate are not 
separated. Moreover, the orbiton-spin interaction has a 
peculiar structure of two fermions interacting with one 
magnon (see below) . To explore the outcome of this com- 
plicated situation, we proceed as follows. Firstly, we work 
well below Tjv. Then the spin operators may be treated 
by the linearized Holstein-Primakoff transformation: as- 
suming a two sublattice bipartite antifcrromagnetic or- 
der, with moments along the magnetic ^-direction, we 
write 



j.mm c ,m 

L ij — °mm' b ij ; 



where = t and differs from zero only when 

(ij) is along the z axis and m = d yz or d zi 
(ij) is along the y axis and m = d xy or d y 
(ij) is along the x axis and m — d zx or d x . 



(2) 



(3) 



The resulting exchange Hamiltonian, for a bond along 
the z-direction, reads 



— Si • Sj 

- 



1 

+ 4 



B,Bj 



B)B\ 



(9) 



Here, Bj is a boson operator which creates a magnon 
excitation (related to Six + iSi V ). Secondly, the KK 
exchange Hamiltonian, together with Eq. (9), will be 
treated as a perturbation, acting on the zeroth order 
Hamiltonian. The latter consists of the free orbiton 
Hamiltonian, and the free spin Hamiltonian. For the 
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first, we impose a constraint ensuring that each site is 
occupied by a single charge carrier. This constraint is 
applied globally. For the latter, we introduce an infinitcs- 
imally small staggered magnetic field in the magnetic z- 
direction, hi = {—l) l h. Hence, 



T~io — Ho + 7~tsi 



(10) 



with 



Ho = -X^(a\a i + b\b i + c\c^, (11) 

i 

Hs = Yl h i(l- B l B i)- (12) 

i 

Note that the chemical potential A includes the constant 
term of the original KK Hamiltonian (7) [cf. Eq. (8)]. 
In the next section we use these definitions in order to 
construct a thermodynamic perturbation theory in terms 
of the Matsubara Green functions. 



III. PERTURBATION EXPANSION FOR THE 
ORBITON SPECTRUM 

The orbiton spectrum is determined by the self-energy 
of the orbiton Green function. Here we calculate the or- 
biton Green function, which is diagonal in the orbiton 
operators, that is, 



G^(r) = (T T a t ,(r)al(0)) 



(13) 



where i and j denote lattice sites, T T is the time-ordering 
operator, and a M = a, 6, or a This Green function satis- 
fies the Dyson equation 

where go is the 'bare' orbiton Green function, 

9o(e m ) = (*£m + A) -1 , (15) 

with e m being the Matsubara frequency. As we show 
below, the self-consistent calculation of the orbiton self- 
energy, M£ v ', requires knowledge of the magnon Green 
function, which may be represented by the matrix 



A,(r) 



(T T Bi(T)B}(0)) (T^WB^O)) 
(TrB^BUo)) (T T Bt (r ) Bj .(0)> 



(16) 



Unlike the standard case of a linear fermion-boson in- 
teraction, the KK Hamiltonian (7) contains six-tail in- 
teraction vertices, of four types (see Fig. 1). These are 
denoted V\, V2, V3, and V4 and are given by 



Vi = -JBjB.al^a^al^, 
V 3 = -JB\B\a\ a tll a\ v a iv , 



V2 has the same form as V± with i replaced by £, and 
V4 — V3 . Hence the KK Hamiltonian involves a fermion- 
magnon interaction, which is quadratic in the magnon 
operators and quartic in the fermion operators. This fea- 
ture is the basis for our results in the following. 
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FIG. 1. The four magnon-orbiton vertices. The wavy lines 
denote the interaction, full (dotted) lines denote the orbiton 
(magnon) operators. 



A. Lowest-order mean-field approximation 



To first order, the perturbation expansion corresponds 
to the mean-field approximation for the self-energy of 
both the orbitons and the magnons. Then, these self- 
energies can be obtained by decoupling the product of 
magnon operators and orbiton operators appearing, e.g, 
in Eq. (4). This procedure is carried out as follows. Ig- 
noring a constant, one first writes 



X (v) Y (v) 



J2[(Y)X {ij) + (X)Y^], (18) 



<«). 



(17) 



where the orbiton operators Y and the magnon opera- 
tors X are given in Eqs. (6) and (9), respectively. The 
notation () stands for the thermodynamic average of the 
relevant operators. A similar decoupling is written down 
for H x and TL V . When these three approximate Hamil- 
tonians are summed together, the resulting Hamiltonian 
separates into two terms, Tii and TL2, the first pertaining 
to the magnons, and the second to the orbitons. 

Let us first consider Hi and the magnon dynamics it 
implies. The orbiton operator product contained in (V) is 
further decoupled, by introducing the thermal averages, 
which will be determined self-consistently (see below), 

n M = ( a ln a i^)^ witn n = n a = n b = n c = 1/3, 

A m = (<4 with A ee A a = A fc = A c , i ^ j. (19) 

In other words, we assume an equal average occupation 
of all three orbitals (i.e., no orbital long range order), and 
consider only the possible "ordering" of the off-diagonal 
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orbital parameter, A, which represents the delocalization 
of the orbital degrees of freedom. Then Hi takes the form 

H i = \ E { B l B * + B ) B i + BiB i + B l B }) > (2°) 

(v) 



where 



with the coupling constant 

J = 2J(n 2 -2A 2 ). 



(21) 



Thus, Hi reduces to the usual spin-wave antiferromag- 
netic Hamiltonian, with a renormalized exchange coef- 
ficient. As usual, this Hamiltonian is diagonalized in 
Fourier space, utilizing the Holstein-Primakoff rotation: 



Cq/3q + Sq(3^_ Q , 



where 



C q = coshWq, 5q = sinhW q 



such that 



■tanh2W, 



q - - (cos q x + cos q y + cos q z ) 



(22) 
(23) 

$q. (24) 



Hence, the magnon Green function takes the form 

£>q(^n) = {iUn - W q ) _1 , 

with the magnon frequency 

Wq = 3 JyJ\ 



$2 

q - 



(25) 



(26) 



Within the mean field approximation of Eq. (18), the 
energy scale of the orbiton Hamiltonian H2 is set by (X), 
which can be now obtained using the magnon dispersion 
(26), 



(X) 



3 
2N 



n- 



+ (2N^ + l)yf 



1 



$2 

q 



(27) 



with the magnon occupation numbers iV q = (exp(/3w q ) — 
l) \ Since the operator still involves products of 
four operators, one needs to apply further contractions. 
It is then convenient to calculate the orbiton self-energy 
using the diagrams shown in Fig. 2. From these diagrams 
one finds the first order self-energy 



M M (k)=4J(X)[n-$CA], 
with the form factors 

— cos k y + cos k z , 



(28) 



(29) 



and analogous expressions for <!> b (y, z — ► z,x) and $ c 
(y, z — > x,y). It therefore follows that the orbiton Green 
function, in this approximation, becomes 



GT(e) 



(30) 



E£ = 4 J(X) [n - $£A] - A, fj, = a, b, c. 



— 1^"-"""^""""^ /vwvwJy' — 



/wv-^-w< 



(31) 
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FIG. 2. The first-order diagrams contributing to the 
self-energy of the orbiton Green function. 



It is thus seen that a finite value of A implies a fi- 
nite dispersion of the orbiton modes, with an orbiton 
band width 16J(X)A. Moreover, this band-width is 
temperature-dependent, due to the temperature depen- 
dence of the magnon occupations, contained in (X), Eq. 
(27), and due to the temperature dependence of A it- 
self. For a given wave-vector k, the three energies 
are generally different, so that the degeneracy of the t2 g 
band is removed. Indeed, were we to terminate the calcu- 
lation at this level, keeping only first-order contributions, 
then the orbiton averages, Eqs. (19), would have to be 
found self-consistently using the Green function (30). It 
is very illuminating to examine those self-consistency re- 
quirements. Equation (30) yields 



(32) 



Then the self-consistency equations read 

1 = ^EEtf. A ( r ) = ^EE*iX (33) 



fi k 



k 



where the first equation fixes the chemical potential A(T), 
so that the three- fold-degenerate orbiton band is 1/3 
filled, and the second determines A. 

Evidently, Eqs. (33) are satisfied by setting A = 0, 
(since then the Fermi function is independent of k, and 
Sk ^k vanishes)- Let us now examine at which temper- 
ature A begins to deviate from zero. To this end we ex- 
pand Eqs. (33) in powers of A. The first self-consistency 
equation yields (at A = 0) 

4J(X)(T)n A(T) = ^J{X){T) - A(T) = ±ln2, (34) 

while the second allows a non-zero A only for tempera- 
tures lower than the solution of 



(35) 
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Using here the explicit expression for (X), Eq. (27), and 
the results Eqs. (21) and (26) for the magnon frequency, 
we find 



2(iJ 



-1 



2(3 J N 



JrE^ coth^. 



(36) 



FIG. 3. Second-order diagrams contributing to the 
self-energy of the orbiton Green function. 

A lengthy calculation of the contribution of the two 
diagrams displayed in Fig. 3 yields (after performing the 
analytical continuation for the Matsubara frequencies) 



This is an implicit equation for the temperature at which 
A begins to deviate from zero. Examination of this re- 
sult reveals that this occurs at ksT of order J. A sim- 
ilar order of magnitude is found from the spin-wave ap- 
proximation for the Neel temperature (above which the 
staggered magnetization (| — B[Bi) vanishes). Both of 
these estimated temperatures are probably above the real 
Neel temperature, which is determined by small symme- 
try breaking terms which are not contained in the simple 
KK model. Hence, in the whole range below the real Tjv 
we expect to have a finite A(T), and we do not expect 
any phase transition at which A vanishes while the spins 
are still ordered. Even if this statement is not valid for all 
T < Tjy, we still expect it to hold for sufficiently low tem- 
peratures. In the next section we examine the dynamical 
corrections to the orbiton spectrum, brought about by 
the second-order corrections, and examine whether and 
how the above conclusion is modified. 



B. Beyond the mean- field approximation 

The next order contribution to the orbiton self-energy, 
M2 j([1 (k, e), requires expansion to second order in the KK 
Hamiltonian, i.e. to second order in the diagrams listed 
in Fig. 1 . Instead of going through the details of all the 
diagrams which contribute to M 2jAi , we give here the de- 
tails of only two diagrams, which describe orbiton prop- 
agation assisted by two-magnon processes (Fig. 3). The 
first of these two diagrams is generated by the vertices 
V\ and Vi in Fig. 1, whereas the second is generated by 
the vertices V3 and V4. Both diagrams contain a single 
propagating fcrmion line. Each of the other two fermion 
tails (of each vertex) is rolled into a bubble, to give the 
factor n 2 . Other diagrams involve three internal fermion 
lines, but do not change the qualitative features discussed 
below. 



N 



u q,p+q 



+ 



7V q (l + 7V p+q ) + / k %(iV p+q - JV q ) 

1 / Ayvp+g + zLpq + iVp+q + Ag 

2 e -^_ p + Wq + Wp+q 

(1 + JV q )(l + JV p + q ) - /£-p(l + ^P + q + JVq) ' 

e - E£_ p -uq- uj p+ q 



for the first diagram, and 
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Nq(l + N p+ q) + f£_ p (N p+ q Nq) 



q,p+q 



k-p + W q - W P+q 



1 + / ^VVp+q + /£-p(l + ^P+q + Nq) 

2 e -££_ p + Wq + u; p+q 

(1 + iV q )(l + AWq) - /£-p(l + ^P + q + ^q) 
£ - K-p - - Wp+q 



(37) 



(38) 



for the second diagram. Here we have introduced the 
magnon coherence factors, 

Cq,p+q = (CqCp+q + SqS p+ q) 2 , 

Cq,p+q = (SqCp+q + CqSp+qf , (39) 

which, upon using Eq. (23), take the simpler form 



C. 



q,p+q 



f cosh 2 (M/ q + VF p+q ) 
\ sinh 2 (W q + W p+q ) 



(40) 





With the second-order corrections to the orbiton self- 
energy, the orbiton Green function becomes 



GT(e) 



1 



(41) 



where M 2 , M (k, e) = M^(k, e) + Af£J(k, e) + .... We now 
concentrate on the contributions of Eqs. (37) and (38) 
to this Green function. Clearly, these contributions yield 
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a finite, temperature-dependent imaginary part to the 
orbiton dispersion, namely, the orbiton modes acquire a 
temperature-dependent damping, as expected. The new 
dispersion law is obtained from the poles of the Green 
function (41), 

z-E£ = M 2 ^(k,z), (42) 

where z is a complex variable. Examining Eqs. (37) 
and (38) one may draw the following conclusions, (i) 
The real part of the second-order self-energy, i?(k, e) = 
Re M 2 , M (k, e), as function of e, contains three compo- 
nents which come from the two-magnon processes. The 
lowest one arises from the two-magnon emission pro- 
cesses and is dominant in the interval — 2ud < e < 0, 
where uod is the maximal magnon frequency [see Eq. 
(26)]. The second results from the processes in which 
one magnon is emitted and one is absorbed, which oc- 
cur in the interval — u D < e < u D . The third compo- 
nent comes from the two-magnon absorption processes 
and prevails at positive frequencies, < e < 2u>u. It 
follows that i?(k, e) has a significant magnitude in the 
entire range — 2ud < e < 2ujd, and falls down rapidly 
outside this interval. The magnitude of different contri- 
butions to i?(k, e) is regulated by the magnon occupation 
numbers. In particular, the mixed emission-absorption 
processes are frozen out at T — > 0. One may conclude 
that the superposition of the three types of the two- 
magnon processes results in a complicated behavior of 
i?(k, e). This quantity may have several zeros and peaks 
within the above energy interval; the heights of these 
peaks can be estimated as ~ kJ, with a numerical co- 
efficient k < 1. (ii) The imaginary part of the orbiton 
self-energy, which comes from the second order terms, 
T(k, e) = Im M 2 ^(k, e), behaves similarly to the real 
part; expressions (37) and (38) do not show that T is 
essentially smaller than R in the whole energy interval, 
(iii) As a result, the orbitons become ill-defined, heavily 
damped, quasi-particles. Their off-diagonal orbiton pa- 
rameter A, and the related orbiton band-width, should 
be again determined self-consistcntly, using the refined 
orbiton Green function Eq. (41), 

A ( T ) = " E E JL f^ ImG r (43) 

in place of the second of Eqs. (33). 

Without going into the details of the solutions of the 
complicated self-consistent equation for A(T), Eq. (43), 
it is clear that going beyond the simple mean-field calcu- 
lation changes drastically the nature of the orbiton ex- 
citations. Whereas in the first order we have obtained 
a coherent orbiton spectrum with a well-defined wave 
vector (appearing due to the magnetic order), this pic- 
ture changes drastically with the second-order correc- 
tions. The interaction with the spin wave excitations ex- 
tends the effective width of the orbiton band and makes 
it asymmetric relative to its center of gravity. Moreover, 



the orbiton excitations lose their coherence by acquiring 
a considerable life-time (given by the imaginary part of 
the self-energy). Those excitations are thus transformed 
into a structureless "liquid-like" continuum. Clearly, 
more complicated processes involving higher-orders of the 
orbiton-magnon interaction, are not expected to modify 
this picture. 



IV. CONCLUDING REMARKS 



We have proposed in this paper a mechanism capa- 
ble of lifting the degeneracy of the i 2 g orbitals in the 
titanate. In our scenario, the orbital degeneracy is re- 
moved due to the interplay between the (ordered) spin 
and the (disordered) orbital degrees of freedom. We 
have obtained three dispersive orbiton bands, with a 
temperature-dependent band-width and a structureless 
spectrum, that exist together with the coherent spin fluc- 
tuations at low temperatures. We have found that the 
second-order perturbation corrections transform the co- 
herent orbiton spectrum, E£, [see Eq. (31)] obtained 
within the simplest mean-field theory, into an incoherent 
continuum of heavily damped excitations, which may be 
termed an 'orbital liquid'. This, in turn, may explain the 
experimental observation, 10 that the low-temperature 
heat capacity arises from the spin wave excitations alone. 

The incoherent continuum state of the orbital exci- 
tations is dominated at low temperatures by the two- 
magnon emission peak in the orbiton self-energy. The 
effect of this peak is to shift the weight of the orbiton 
spectrum into negative energies, thus making the orbiton 
liquid state thermodynamically preferable as compared 
to the fully-degenerate orbiton state that prevails when 
A = 0, or in the absence of the coupling with the spin 
degrees of freedom. In particular, our process of de- 
generacy lifting may possibly overcome the Jahn- Teller 
mechanism. 13 In the latter, the energy gain is restricted 
by the energy of the Jahn- Teller coupling strength, which 
is definitely smaller than the exchange coupling, which 
determines the energy gain in our case. 

Finally we mention that, within the framework of the 
KK Hamiltonian in cubic symmetry, the orbiton-magnon 
interaction cannot lead to a magnetic anisotropy, that 
is, it cannot generate a gap in the spin-wave spectrum. 
To obtain such anisotropics, one should extend the KK 
Hamiltonian to include the hybridization of the oxygen 
and titanium orbitals, the possible orthorhombic tilt- 
ing of the oxygen octahedra, as well as spin-orbit in- 
teractions. Apparently, one will also have to treat the 
Coulomb interaction in a more refined way, taking into 
account Coulomb-exchange effects. 5 
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